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Abstract 



We consider homogeneous multidimensional continued fraction algorithms, in partic- 
ular a family of maps which was introduced by F. Schweiger. We prove his conjecture 
regarding the existence of an absorbing set for those maps. We also establish that their 
renormalisations are nonergodic which disproves another conjecture due to Schweiger. 
Other homogeneous algorithms are also analysed including ones which are ergodic. 

1 Introduction 

Multidimensional continued fraction algorithms have become a classical topic in the Theory 
of Dynamical Systems. There is a wide literature on the metric theory of those maps and a 
large collection of examples is available in Schweiger |10j . Besides their number theoretical 
motivation they also appear in dynamical systems, in particular within the renormalisation 
theory where Poincare first return maps are considered. A seminal example is the so-called 
Rauzy induction algorithm of interval exchange transformations [9] which is a central tool 
to study the ergodic properties of interval exchanges, e.g. Veech [111 I12j . Nogueira and 
Rudolph |8j and Avila and Forni [TJ. The Rauzy algorithm relates one interval exchange 
to a class of interval exchange maps through a suitable induction process. 

The present work concerns the so-called homogeneous subtractive algorithms. One of 
the difficulties that one encounters in studying their dynamics lays in the fact that they 
admit an infinite invariant measure. Most of the algorithms which have been studied bare 
ergodic properties, less attention has been paid to nonergodic algorithms having absorbing 
set. Homogeneous algorithms having absorbing set are quite common maps though. Here 
we consider a class of maps which are naturally defined and prove that they have absorbing 
set. 

Among the well known examples of homogeneous maps we enumerate the Euclidean 
algorithm, the Brun algorithm and the Selmer algorithm. Quite often when a nontrivial 
absorbing set occurs a suitable renormalisation of the algorithm is considered forcing its 
ergodicity, but very little is known about the dynamical properties of the initial homoge- 
neous algorithm. 

Here our purpose is to study the family of transformations defined as follows: Let 
A n = {x = (xi, X2, ■ ■ ■ , x n ) : < xi < . . . < x n } be the set of those points of W 1 whose 
coordinates are positive and nondecreasing. Let a,b > 1 be integers. Given any point 
x G A a+b , we define 



where ir is a permutation of the coordinates (depending on the point x) which arranges 
an, ... , x a ,x a+ i — x a , . . . , x a+ i, — x a in ascending order. The map T a ^ : A a+b — > A a+b is 
continuous and piecewise linear. 





The family (1.1) is not new and contains several well studied maps given by special 
values of the parameters a and b. The first example, corresponding to a = b = 1, has 
the following property: If x = (xi,X2) G A 2 n Z 2 , then there exists k > 1 such that 
T^ 1 (x) = (0, d), where d > 1 is the greatest common divisor of x\ and x%. That explains 
why it is called the Euclidean algorithm. It is known that for almost every x G A 2 the 
corresponding orbit converges to the origin and the map T\ t \ is ergodic with respect to 
Lebesgue measure [T2"l 17]. 

Another special case corresponds to a > 2 and 6=1 which is called Brun algorithm. 
Later it will be shown that for almost every x G A a+1 the corresponding orbit also converges 
to the origin. In order to study the metric properties of this algorithm, one may project 
the dynamics onto the set B = {x G A a : < x\ < . . . < x a < 1}, which is of finite 
Lebesgue measure. The new map S a : B —> B is defined through the commutative diagram 

T 

^a+l ft' 1 v , ^a+1 
V 

B ► B 

S a ,i 

where p is the projection (xi,x 2 , . . . ,x a+ i) H> (xi/x a+ i, x 2 /x a+ i, . . . , x a /x a+ i). A result 
by Schweiger ([E], p. 50) states that S a ,i is ergodic with respect to Lebesgue measure. 
Nevertheless the question of ergodicity of the original map T 0) x : A a+1 — > A a+1 remains 
open. 

The third type of map corresponds to a = 1 and b > 2 and was considered by Meester 
and Nowicki [5J and later by Kraaikamp and Meester [3]. The map appeared naturally 
in their study of a parametric percolation model on the lattice Z fe+1 . They proposed an 
algorithm for computation of the critical probability of the model. The efficiency of this 
algorithm may be expressed through the absorbing properties of the set 

A = {x G A b+1 : xt + x 2 + . . . + x b+1 < bx b+1 } 

under the dynamics of T = T\ We recall that T is an absorbing set if almost surely 

|J T'\A) = A 5+1 . 

To be more precise, if A is absorbing, then the algorithm gives the critical probability of 
the model in a finite number of steps, for almost every choice of model parameters. The 
case b = 2 was first studied in [5]. In [3J the results were extended to the case b > 2. 

The general case of the map T a ^ : A a+b — > \ a+b was first studied by Schweiger (see 
|10] , chapter 9). He noticed that for every a,b> 1, the corresponding set 

A = {x G A a+b : xi+x 2 + ... + x a+b < bx a+b } 

was T a>b - forward invariant and conjectured that this set should be absorbing under T a ^ b . 
Moreover, in the case b > 2, he defined the smaller invariant set 

D = {x G A a+b : xi+x 2 + ...+ x a+1 < x a+2 } 

and conjectured that it was absorbing as well. 

Conjecture 1 (Schweiger). Let a > 1 and b > 2. Then the corresponding set D is 
absorbing for T ajb . 
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In order to prove this conjecture, we first establish an extension of the Kraaikamp- 
Meester-Nowicki theorem. 



Theorem 1.1. Let a > 1 and b > 2. The set A is absorbing for T a ^. 



Analysing the limiting behavior of the orbits, we prove a stronger claim. 



Theorem 1.2. Let a > 1 and b > 2. Then, for almost every x £ A' 



limT a fc 6 (x) = (0,...,0,x- 2 ,...,x- b ), 



where < x 




< . . . < x 




As an immediate corollary, we get that Conjecture 1 holds. 

Schweiger also attempted to study the ergodic properties of the map T a ^. To this end, 
he used the same projection as in the case of the Brun algorithm. The new transformation 
S a 6 : B — > B is defined through the commutative diagram 



where B = {x G A a+b 1 : < xi < . . . < x a +&-i < 1} and p is the projection 

(X!,X 2 , ■ ■ ■ ,X a+b ) H-> (x 1 /x a+b ,X 2 /x 

a+bi • • • i x a+b-i / x a+b) ■ Let D stand also for the pro- 
jection of the set D C A a+b onto B. This projected set is S^b-forward invariant, just 
as the original D is for T a ^. The transformation S a ^ restricted to the set D admits an 
invariant measure which is absolutely continuous with respect to Lebesgue measure and 
whose density was explicitely calculated by Schweiger ([TO], p. 80). Having this in mind, 
it is natural to ask about the ergodic properties of S a j> with respect to Lebesgue measure 
on D. The following conjecture was formulated. 

Conjecture 2 (Schweiger). Let a > 1 and b > 2. Then the restriction of the map S at b to 
the set D is ergodic with respect to Lebesgue measure on D. 

This conjecture is more subtle. First, we prove the following theorem that confirms 
the conjecture in the case a = 1 and 6 = 2. 

Theorem 1.3. The map Sip : B — )■ B is totally dissipative and ergodic with respect to 
Lebesgue measure. 

On the other hand, from Theorem |1.2| we deduce that for almost every x E D we have 
limjt^oo S% >b (x) = (0, . . .,0,a;^. 2 , . . . ,x ! ^ +b _ 1 ). If b > 3, we may show that f(x) = x~ is 
a nonconstant invariant function for S ajb . This disproves the second Schweiger conjecture 
in this case. 

Theorem 1.4. For b > 3, the map S a ^ is not ergodic with respect to Lebesgue measure. 

The same argument shows that for any a > 1 and b > 2 the nonprojected algorithm 
T a £ ■ A a+b — > A a+b is not ergodic with respect to Lebesgue measure. The question of the 
ergodicity of the projected transformation S a ^ for a > 2 seems more difficult and is not 
answered here. 
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The paper is organized as follows. The next section regards a transformation of 
which is closely related to the map T± ; 2 of our family. We prove the first Schweiger con- 
jecture in this special case. The Section 3 is devoted to some general remarks concerning 
the family of transformations T a ^. Some notations used throughout the text are also in- 
troduced. In Section 4, which is highly inspired by the work of [3], we characterize the 
invariant sets A and D through the limiting behaviour of the orbits. In order to show that 
the set A is absorbing, we have to study the dynamics of T a ^ on its complement which is 
done in Section 5. In particular, we define some smaller set contained in the complement 
of A and define a kind of first return map induced by T a ^ on this set. This is done despite 
the infinite measure of the sets involved. A family of partitions of 0, connected to the 
iterations of the first return map, is defined. In Section 6, we project the dynamics onto 
a finite Lebesgue measure set of lower dimension which enables us to study metric prop- 
erties of the partitions mentioned before. Those properties are used in Section 7 to prove 
Theorems |1.1| - 11.4[ In the last section we define other families of subtractive algorithms 
and discuss the possibility of extending our results to those new transformations. 



After we concluded the present paper we have been informed that an independent 
article [2J with similar results (Theorem 1.2 and 1.4) is currently accepted for publication 
by Israel Journal of Mathematics. However, we would like to mention that our paper 
presents a different approach to the subject and places it in a larger perspective. In 
particular it allows us to prove Theorem |1.3| 



2 An example 

The present section deals with a particular version of the map and is independent of 
the remaining of the paper. It illustrates our geometric approach to the proof of Conjecture 
1, which will further be extended to the general case of the map T a ^. 
Let T : R 3 + -> R% be given by 

{(xi,x 2 - 21,23 - 21) if min{2i, x 2 , 23} = x\ 
(21 - 22,22,23 - 22) if min{2i, 22,23} = 2 2 
(21 - 23,22 - 23,23) if min{2i, 22,23} = 23 

which is well defined for almost every x G We will show that the set 

A = {(21,22,23) G M.+ : 21 + x 2 + 23 < 2max{xi,2 2 ,2 3 }} 

is absorbing for T which will clearly imply the Conjecture 1 for the map 7i 2- An argument 
similar to the one presented below, due to J.-C. Yoccoz, can be found in [3J. 

Let us introduce the following notation. Let F = (/i,/ 2 ,/3) be a basis of M 3 . We 
define the corresponding positive simplicial cone 

M +(- ?r ) = {xeR%:x = 21/1 + 22/2 + 23/3, 21, 22, 23 > 0} 

and the analogue of the set A relative to this new basis 

A(F) = {2 = 21/1 + 22/2 + 23/3 : (21,22,23) G A}. 

The set A{F) is equal to the union of three cones M+(/i,/i + /2,/i + h), K+(/2,/i + 
f-2, h + h) and R% (/ 3 , f x + / 3 , h + h\ 
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A simplicial cone R+(/i, /2> /k) C may be visualized as its intersection with the 
2-dimensional simplex A = {x G : x\ + X2 + x 3 = 1}. The cone is then represented by 
the triangle whose vertices are the radial projections of vectors /1, f 2 and /3 onto A. For 
convenience we choose to label these vertices with the original coordinates of the vectors 
fi rather then their projections. 

In particular, the original set A may be seen as A(ei,e2,e 3 ), where (ei, 62,63) stands 
for the canonical basis of M 3 . We get that A is a union of three simplicial cones 

A = A(ei,e 2 ,e 3 ) = M 3 (ei,ei + e 2 ,ei + e3)UM 3 (e2,ei + e2,e2 + e3)UlR+(e3,ei+e3,e2 + e3). 

Its image in the simplex A is the union of three shaded triangles in Figure 1. 

Let now x = (xi, X2, X3) represent the vector x in the canonical basis of M 3 . Consider 
the cone M+(ei + e2 + e 3 , e%, e^) corresponding to the set of vectors such that x\ < X2, x 3 . 
We remark that 

x = x\e\ + x 2 e 2 + x 3 e 3 = x\{e\ + e 2 + e 3 ) + (x 2 - x\)e 2 + (x 3 - xi)e 3 . 

This relation means that the map T may be seen as a convenient basis change in M 3 and 
implies 

T~ 1 (A) n R + (e a + e 2 + e 3 , e 2 , e 3 ) = A(e x + e 2 + e 3j e 2 ,e 3 ). (2.1) 

The projected version of this relation is depicted in Figure 2. The analogous relations hold 
in M 3 (ei + e2 + e 3 , e\, e 3 ) and M 3 (ei + e 2 + e 3 , ei, e 2 ) where x 2 < xi, x 3 and X3 < xi, x 2 
respectively. Putting them together yields 

T~ l {A) = A U M^(2ei + e 2 + e 3 , ej + 2e 2 + e 3 , ei + e 2 + 2e 3 ). 




(|,0,|) 



(1,0,0) 



(iio) 

Figure 1 



(0,U) 



(0,1,0) 



(0,0,1) 




Figure 2 



(0,|,|) 



(0,1,0) 



In order to study the second inverse image of T, we take the cone Ki(ei + e 2 + e3, 62, 63) 
and divide it into three sub-cones M5j_(ei+2e2+2e3, e 2 , 63), ^^(61+262+263, ei+e2+e3, 63) 
and M+(ei + 2e2 + 2e3,ei + e2 + 63,62), which correspond to the choice of the smallest 
coordinate with respect to the basis (ei + e2 + 63, e2, 63). We get the following counterpart 



of the relation (2.1 ) 



T 2 (A) n M 3 (ei + 2e 2 + 2e 3 , e 2 , e 3 ) = A{ei + 2e 2 + 2e 3 , e 2 , e 3 ) 

and the other two corresponding to the remaining two sub-cones. 

Arguing by induction, one may show that in the kth step we obtain the decomposition 
of into 3 fc simplicial cones, with disjoint non empty interiors. If M^(/i, f 2 , / 3 ) is a cone 
of the kth decomposition, then 

T-\A) n M 3 (/ l5 f 2 , f 3 ) = A(h, f 2 , f 3 ). 
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In the next step, we decompose R^(/i, f 2 , f 3 ) into three sub-cones corresponding to three 
basis (/1 + h + h, fuh), (fi +f 2 + f 3 , f2, f 3 ) and (/1 + f 2 + f 3 , fi, f 3 ). 

We may also describe this construction by putting the stress on the complement of 
the inverse images of A. For every k > 0, the complement in of the set T~ k (A) is 
composed of 3 k simplicial cones of disjoint interiors. If K^_(/x, f 2 , f 3 ) is one of them, then 

T^ k+l \A) n R 3 + (h, f 2 , f 3 ) = Rl(fr + f 2 , f 2 + h, h + /s)- (2-2) 



(- 1 ) fi 1 M (0 1 i) 



/l 1 In /ll IN 

V 3 ' 6 ' 2> Vk'3'2' 1 




(- 1 M 

V 6 ' 2 ' 3 ' 




Figure 3 



Figure 3 depicts the first three steps of this construction projected onto A and restricted 
to the "middle" triangle of vertices ^,0), (0, |, |) and (^,0, ^). 

In order to show that the set A is absorbing, we have to show that, for every step 
k > 0, the cone + f 2 , f 2 + fs, fi + f 3 ) corresponds to a large part of the volume of 

the cone f 2 , f 3 ) in (2.2). Since both cones are of infinite measure, we compare the 



areas of their images in A. Let &-f 1 f 2 f 3 be the triangle in A corresponding to the cone 
Hi f 3 )- It may be shown (see [6], Lemma 3.2) that 



Area(A flj2 j 3 ) 



V3 \det(M fl)h)h )\ 



2 H/1II1II/2IUII/3II1' 

where Mj x j 2 j 3 is the non negative matrix whose columns are given by the coordinates of 
the vectors /1, f 2 and f^. 

Let k > and Af lt f 2 j 3 be one of the triangles in the complement of T~ k (A). By 
induction we show that any of the mutual ratios of the numbers ||/i||i, H/2II1, II /3 111 does 
not exceed k + 1. It is clear in the case of the unique triangle in the complement of A 
(k = 0). Assume the assertion is true for the triangle A/ x j 2 j 3 of the complement of 



T k {A). This implies the assertion for the three triangles A/j fi + f 



2,/l+/ 3 ) A /2,/l+/2,/2+/3> 



^/3,/i+/3,/2+/3 °f the complement of T ^ + '(A). Indeed, for Af 1 j 1+ f 2 j 1 +f 3 we get 

II/1 + / 2 ||i < IIA + /2II1 =1 { II/2II1 <k | 2 



II/: 
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II/1 + /3II1- H/llll 
Applying this estimation we get 

Area(A h+hj2+f3jl+f3 ) __ 2||/i||i||/ 2 ||i||/3||i 



Area(A flj2j3 ) 



> 



II/1 + /2II1H/2 + /3II1II/1 + /3II1 " 2(k + 2)- 



This together with (2.2) implies that 



Area(A \ T^ k+1 \A)) < (l - ^^y) Area(A \ T~ k (A)) 
for every k > 0. As the series J2 k>Q 2(fc+2) diverges, the set A is absorbing for the map T. 
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3 Some notations and remarks 



Throughout the remaining part of this paper the following notations are used. All elements 
of W 1 are given in row notation. The canonical basis of R n is denoted by (ei, e2, . . . , e n ) 
and all vectors of R n are always expressed with respect to this basis. Accordingly, the 
expression x = (x%, X2, ■ ■ ■ , x n ) represents the vector x in the canonical basis. Where 
matrix multiplications are used, the symbol x T stands for the transpose of x, i.e. the 
column vector corresponding to x. We also define <r(x) to be the sum of the coordinates 
of a vector x, that is if x G W 1 then a{x) = x\ + . . . + x n . The actual value of n will be 
clear from the context. 

While studying a transformation T : W l — > 1" we will often write x^ instead of T k (x), 

(k) 

for x G W 1 and k > 0. Accordingly, for any i = 1, . . . , n, the symbol x\ will stand for the 
ith coordinate of the point x^ k \ 

When speaking about a measure, we always mean Lebesgue measure in the corre- 
sponding M n space. Thus all "almost surely" and "almost every" statements refer to the 
corresponding Lebesgue measure. 

Let a, b > 1 be integers. Let x G A a+b and define y G by 

V — (*El> • • • j X a > ■Eq+I -^d) • • ■ 1 X a -{-b %a) • 

For almost every x there exists a unique permutation 7r x of the set {1,2,... ,0 + 6} such 
that {yn x (i), yn x (2), ■■•> U-K X {a+b)) € A - a+b - The onl y ambiguity appears when two or more 
coordinates of y are equal. In such a case, we define tt x to be the only permutation which 
satisfies 7r(i) < 7r(j) anytime j/, = y 7 - an d z < j. 



We may thus redefine the map (1.1) as the transformation T a & : A a+b — > A a+b defined 
by 

T a ,b(x) = 1/^(2), •••,l/«- a ,(o+6))- 

The map T & is continuous and piecewise linear, where the number of pieces corresponds 
to the number of different permutations ir x involved. Since for every x G A a+f> we have 

xi < x 2 < ■ ■ ■ < x a and x a+ i - x a < x a+2 - x a < . . . < x a+b - x a , 

every permutation ir x satisfies vr z (l) < . . . < ir x (a) and ir x (a + 1) < . . . < ir x {a + b). The 
number of permutations is thus equal to the number of different shuffles of a deck of a 
cards with a deck of b cards, namely ( a ^ b )- 

Let Tl a ,b be the set of all these permutations. For every tt G H a ,b we define the cylinder 
A n = {x G A a+b : tt x = tt}. The cylinders A,r define a partition of the set A a+b into ( a + 6 ) 
regions, with nonempty and disjoint interiors, which intersect along null Lebesgue measure 
sets. We also define a family {L w ,ir G H a ,b} of m atrices in GL(a + &,Z) such that T a ^ 



restricted to A,r is given by x h-> L w x t . From (3.1) we easily deduce that the line vectors 
of the matrix L n are 

61; 62) ■ • • j 6a— 1) 6a > 6a+l 6 a , . . . , 6 Q -_)_ft 6 a , (3-2) 

rearranged in the order given by the corresponding permutation tt. To be more precise, 
if Lj stands for the jth line of the matrix L n , we have L^u) = &i for i = 1, . . . , a and 
= ej — e a for i = a + 1, . . . , a + b. 
The map T a ^ is not globally injective. The next lemma comes from [10J. 

Lemma 3.1 ([TO], p. 69). Every cylinder A n is full, that is T a ^{A v ) = A a+b . 
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We may also consider the family of inverse matrices {M^, tt G n a ,b}i where M n = L^ 1 , 
corresponding to all inverse branches : A a+b — > A n of the transformation T a ^. For every 
permutation tt the corresponding matrix M n is a non-negative element of GL{a + b, Z). 
From (3.2) we get that the column vectors of M v are 

e ii ■ ■ • i e a-i) i e a + e a +i + . . . + e a+b ) T , ej +1 , . . . , e^ +b , (3-3) 

arranged once again in the order given by tt. If Cj stands for the j-th column of M n , we 
have C w (q = ej for % a and C^a) = ( e a + e a +i + ■■■ + e a+b ) T '. 

We are interested in the absorbing properties of the following sets 

A = {x G A a+b : o(x) < bx a+b } and D = {x G A a+b : x\ + x 2 + . . . + x a+ i < x a+2 }, 

the former being defined for all a, 6 > 1 and the latter for b > 2. In the case 6=1, the set 
A is a null measure set. For b > 2, both sets ^4 and L> are infinite measure sets. Moreover 
D C A with .D = A if and only if b = 2. We recall the following result due to Schweiger: 

Lemma 3.2 ([ID], p. 72 and 78). The sets A and D are forward invariant for T a ^, i.e. 
T atb {A) = A and T afi (D) = D. 

Our goal is to show that they are actually absorbing for T a b . 



4 Limiting behavior of orbits 

It is a trivial remark that the orbit of any point x G A a+b under the action of T a ^ is 
convergent in A a+b . Indeed, for any i = 1, . . . , a + b fixed, the corresponding coordinate 
sequence (x[ k ^)k>o, is positive and nonincreasing. We will show that the set A attracts 
almost every point x 6 A a+b for which x^ do not converge to the origin as k — > oo. 
Equivalently, almost every point whose entire orbit stays outside of A is characterized by 
the condition x^ — > (0, . . . , 0). In the case a = 1 and b > 2 this was shown in [3]. 

Theorem 4.1 ([3], Theorem 1 and 2). Let b > 2. For almost every x G A b+l the sequence 
TJ*Jx), k > 1, does not converge to the origin and the set A is absorbing under T\^. 

We begin with results that illustrate the dynamical behavior of homogeneous algo- 
rithms. 

Lemma 4.2. Let a > 1 and b > 1. For almost every x G A a+b we have 

lim T a u[x) = (xf 3 , . . . , x^ +b ), 

k— >+oo 

with X f = xf = ... = x™ +1 = 0. 

Proof. Let x G A a+b . For every i = 1, . . . , a + b, the sequence (xj^)fe>o> is convergent to 
some We obviously have < x^° < x?f < . . . < x^ +b since at every iteration the 
coordinates are reordered to be nondecr easing. 

Suppose that the coordinates of the point x are rationally independent. This implies 

(k) 

that x- > for all i = 1, . . . , a + b and k > 0. For every A; > 1 there exists m > 1 such 
that x^\ — mx^ < Xa \ If m is the smallest one for which this inequality is satisfied, we 
get x^+™^ = and the two limits x^ and x^ +1 must be equal. 
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Since x$ and x^h converge to the same limit, for every e > there exists k > 1 such 

that < x^h — Xa < e. This implies that x[ k+1 ^ = min{xi > x a+i ~ x a} — e - We get 
xf = 0. 

Suppose that x?° > for some 2 < i < a + 1 and let j = min{2 < i < a + 1 : x?° > 0}. 
Since xj } _ l = 0, there exists some k > 1 such that x^.j^ < x°°/2 and ~~ ^ — x f 
This implies £j fc+1 ' > = max{xj_ 1 , a^+i — si } < Xj° /2, which is in contradiction with the 
fact that the sequence Xj is nonincreasing. We get x°° = for every i = 1, . . . ,a + l. □ 

Another useful result states the following. 

Lemma 4.3. Let a > 1 and 6 > 1. Let x G A a+6 satisfy x^ — > as — > +oo. Taen 
(x(x) = b^2 k>0 x^ and thus the sum involved is finite for almost every x. 

Proof. For any m > 1 we clearly have a{x^ m ^) = a(x^ m ~ 1 ^) — 6xi m ^ which gives 

m— 1 
fc=0 

Since cr(x ( - m * ) ) goes to zero as m goes to infinity, we get the claim. □ 

The last two lemmas applied to the Brun algorithm T aj i : A a+1 — > A a+1 give 

Corollary 4.4. Let a > 2, b = 1 and write T^(x) = x^ k \ Then for almost every x G A a+1 
we get 

oo 
k=l 

We will now describe the generic limiting behavior in the set A. 



lim x^ = and a(x) = xi k 



Lemma 4.5. Let a > 1 and 6 > 2. TTien lim/^oo T k b (x) ^ (0, ...,0) for almost every 
x £ A. 

Proof. First, remark that 

T (*+i) _ max ; T W J*) _ T (fc)\ > J k ) _ T (k) 

for every k > 0. This implies that a^+fi > x a+ ^ — Yli=o Now suppose that x^ — > 

we get cr(x) = b 5^fc>o • If # € -A, we also have x a+ fc > 



4.3 



as A; — )• oo. From Lemma 

a(x)/b = Y^k>o x " ■ Since x^ fc ^ — >■ 0, we must have 

fc-1 

0<x a+6 - J>« <i^->0. 

8=0 

This implies x a+ ;, = o~(x)/b which is a null measure condition. □ 
The proof of the next proposition is a straightforward generalization of an argument 

of m- 

Proposition 4.6. Let a > 1 and b > 2. Assume that the coordinates of x G A a+b are 
rationally independent and suppose that x^ -fr (0, ... ,0) as k — > oo. Taen i/iere exists 
> 1 snc/i taat x W G A. 
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Proof. Let x be such that lim^oo = (0, . . . , 0, x^ +2 , ■ ■ ■ , x^ +b ) ^ (0, . . . , 0) and let 
m = min{z : ^ 0}. Obviously 2 < m < b. Since all the coordinates to the left of 
Xa+m go to zero as k — > oo, there exists k > 1 such that 

(fc) (fc) (fc) 

X 1 j-Xg + ■ . . + X a+m _ 1 (fc) 

m - 1 " Xa+m ' 

Therefore we may write 

T (*) + T (*)+ +T ( fc ) +T ( fc ) T ( fc ) , r W , + r « _ -i _(*) 

m m — 1 m m ~ 

(k ) f m-1 1 \ _ (fc ) 

— x a+m I m m I ~ X a+m — X a+m+l- 

Finally we get 

(fc) (fc) (fc) 

x l + ^2 "+" • • • T X a+fe (fc) 



which means S A □ 



As far as the set D is concerned, we may strengthen the conclusion of Lemma 4.5 
Lemma 4.7. Let a > 1 and b > 2. For almost every x £ D we have 

lim T a fe fe (x) = (0,...,0,x~ 2 ,...,x^), 

fc— >oo 

where x?° > 0, j = a + 2, . . . , a + 6. 

Proof. The property x a+ 2 > xi + . . . + x a+ i implies that for any x E D we have 

^(x) = (T 0) i(»i, . . . ,X a +i),X a+2 - X a , . . . ,X a+b - X a ). 

Since T(D) C F>, iterating the previous formula gives 

fe-i fc-i 



r£ 6 (s) = (r a fc 1 (x 1 ,...,x a+ i),x a+2 -^x«,...,x a+fe -^x«), (4.1) 
for any x £ D and fc > 1. Let now x £ D satisfy T^ +1 (xi, . . . , x a +i) — > and x a+ 2 > 



x\ + . . . + x a +\. Corollary 4.4 implies that this is true for almost every x G D and also 
that x°° = Xj — a(x±, . . . , x a +i) > 0, j = a + 2, . . . , a + b. □ 

5 Dynamics outside A 

In order to show that A is an absorbing set we have to understand the dynamics of T a ^ 
on its complement. Define 

cA = A a+b \A = {x£ A a+b : a(x) > bx a+b }. 

We begin with the following remark. 

Lemma 5.1. Let x 6 cA and suppose that x a +b > 2x a . Then T a fi(x) 6 cA 
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Proof. Let x G A a+b and write y = T a>b (x). We have y\ = min{a;i, » a +i — ^a} and y a+ & = 
max{x a , x a+ b — Xa}. If we suppose that x £ cA and x a+ b > 2x a , then y a +b = ^a+fe — ^a- 
We get 

by a +b = bx a+b - bx a < a(x) - bx a = a(y), 
which means y G cA. □ 

Having this in mind, we define the following subset of cA 

Q = {x £ cA : 2x a > x a+b ). 

The previous lemma implies that any orbit starting in cA and eventually attracted by A 
has to visit 0. First, we will show that almost every orbit starting in cA visits (which 
does not mean it is attracted by A though). 

Lemma 5.2. For almost every x G cA there exists k > such that T^ b (x) G 0. 



Proof. Let x G cA be a point whose orbit never leaves cA. We know (Proposition 4.6) 
tha t in this case x^ converges to the origin almost surely, which also means (Lemma 



4.3) that o~{x) = b^2 k>Q x^ ■ If the orbit of x never visits 0, for every k > we have 



x a+b^ = x a+b ~ x< a \ which implies x a+ b = Ylk>o x ^ = a ( x )/b- This contradicts the 
definition of cA. □ 

For convenience we set 

T = cA^T~ 1 {A) C 0. 

To show that A is absorbing for T a j,, it is enough to show that the set V absorbs almost 
every orbit starting in cA. To this end, we consider a kind of a first return map P : — > 0, 
conditioned on the set T: 

{x if x G r 

T a>b {x) if x T and T a>b (x) G 
Tl b {x) if x £ T, T a fc 6 (x) G and T^x) £ 0, i = 1, . . . , k - 1. 

This map is well defined for almost every point in 0. 

The map P is piecewise linear just as T a ^ is, but now the number of cylinders is infinite, 
since the first return time to is not bounded. Let us define them explicitly. For every 
k > 1 and every set of permutations 7ri, . . . , tt^ G H a ,b, we put 

©Hi,...,** = {x G \ r : P(x) = T^ b (x) and T^ b (x) = L Vi ■ ■ ■ L ni x T for i = 1, . . . , k}. 

Some of these cylinders are empty. If a given cylinder is nonempty, then it corresponds 
to the set of points of 0, whose orbits under T a ^ visit the sequence of cylinders A^, 
i = 1, . . . , k, before returning to the set 0. On such a cylinder, P is given by the matrix 
£7ri,...,7r fc = L nk • • • L 7Tl which is a product of a finite number of matrices L^, corresponding 
to the cylinders of T a ^ visited by the orbit. We may thus decompose as follows 

oo 

= ruj |J (-), (5.i) 

fc=i (*i,...,**)en; j 6 

The cylinders have the following property. 

Proposition 5.3. For every choice tt%, . . . , tt/. G H a ,b the corresponding cylinder & ni „,^ k is 
either empty or full for the map P, i.e. ^(©ti-!,...,^) = whenever 0^,...,^ is nonempty. 
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Proof. The set G is defined by the inequality x a > x a+ b — x a , where both terms are the 
coordinates of T a ^{x). Since every cylinder A^ corresponds to some order of coordinates 
of T a ^(x), for every it G U at b we have either A^ n cA C or A„- n cA C cA \ 0. We have 
the following: 

P(e*i,...,0 = T aA x e \ r : x ^ gA 12 ,., .,x^ g A, fcI x w g ©}, (5.2) 

where all the cylinders involved are restricted to cA. It is obvious that our cylinder is 
empty, unless A vri n cA C and A Wi n c^4 C c0 for i = 2, . . . , k. 

Suppose that our cylinder is nonempty. Recall that the cylinders A^ corresponding to 
the map T a ^ are full. Since T a ^{A) C A and T a ^{T) C A, we must have T a fi(A. n r\(cA\T)) = 
cA for every ir G H a ,b- This, together with (5.2), clearly implies that P(07n,...,7r fe ) = 0. □ 

Let us slightly alter the notations in order to introduce the sequence of refinements of 
the partition (5.1). Since this partition corresponds to the first iterate of the map P (or 
the first return time of T & to 0), we will call its cylinders , where v is some finite 
sequence of permutations from the set n a b. We may rewrite (5.1) as 



© = ruU U 

k=i v&n k , 

a.b 

Every cylinder 0^ may be partitioned further, considering the second iteration of P. 
For every k > 1 and u' G 11^ b we define 

0<L = {x G 0« : P(x) G q% 



which generates the following refinement of the partition (5.1): 

oo 

0=p- i (r )u u u u <e- 

fe,m=i ueu k a /'en™ 6 

It is not hard to show that the cylinders of this new, finer partition are either empty or 
full for the map P 2 . 

Following the same pattern, for every n > 1 we may consider the nth iteration of the 
map P and define the corresponding nth partition of 0: 

oo 

= P -(n-i) (r)u |j |j ... (j eSL^, (5.3) 

fcx,*2,...,fc„=i I/l6n *i j u n en h a " b 



©S...^ n = {x G eg) : P(x) G eg), . . .,P n -\x) G ©«}. 



where 



The following result may be deduced by induction from Proposition |5.3| and the definition 
of the nth partition. 

Proposition 5.4. For every k\ , . . . , k n > 1 and every vi G II A , i = 1, . . . , n, the corre- 
sponding cylinder ©i™!...*^ is either empty or full for the map P n , i.e. F n (©S...^J = 
whenever ©[/"!.. .*„ n is nonempty. 
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6 Projected dynamics 



In this section, we will show that some constant proportion of points from every cylinder 
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(n) 



falls into T under the action of P n . To make this statement meaningful, we 



project our dynamics onto the simplex A = {x G A a+b 
(a + b — l)-dimensional Lebesgue measure. Let p : A a+b - 
x i-> x/a(x) and consider the new transformation T a ^ : A 
diagram commute: 



cr(x) = 1}, which is of finite 
A be the projection given by 
► A which makes the following 



P 

A 



L a,b 



> A a+b 
P 

> A 



l a,b 



We may define the projected counterparts of the sets A and D which are given by 

A = {x G A : 1 > bx a+ t,} and D = {x G A : x\ + X2 + . . . + x a +i < x a +2}- 

It is easy to see that both sets are T a) b-forward invariant. We will also consider the 
projected set 

6 = {x G A : 1 > bx a+b , 2x a > x a+b } (6.1) 

with its subset T, the corresponding first return map P : Q — > Q and the underlying 
partitions into cylinders G>i"i...^ n . 

Proposition 6.1. There exists a constant a G (0, 1) such that for every n > 1 and every 
projected cylinder 6i™i...«, n of the nth partition of @, we have 

~,eb(x F B, ( "l ... : P n (x) G f) 

> a, 



where Leb stands for {a + b — 1)- dimensional Lebesgue measure in A. 



(6.2) 



Proof. Let G*!™!...*^ be a cylinder. As the initial cylinder 0i,i*...*i/„ is full for the map P n , 
we easily get that P n (0 



(n) 



(n) 



l>l*...*V n 



' (n) 

0. Moreover, on 0^4...^ the map P n is defined by 
Lx T 



P n (x) 



a(Lx T ) ' 



(6.3) 



where L is the finite product of matrices L Wi corresponding to the cylinders of the original 
transformation T a ^ (the projection map p : A a+b — > A preserves the order of coordinates 



at any point of A a+6 ). Let M = L~ l and denote by P^ n : 
inverse branch of P n . Once again we have 







Mx T 



We get 



Leb(Q 



Vl*...*V>, 



a(Mx T ) 

■«)=/-„ ^=/- - 



the corresponding 
(6.4) 



/_ J(P M n )dx 
Je 
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and 



Leb(x e 6 



(n) 

Vl*...~kl/ n 



: P n (x) G f ) 



dx 



J(PM n )dx, 



where J(P^ n ) stands for the Jacobian of the transformation defined by (6.4) and dx is 
(a + b — l)-dimensional Lebesgue volume element on A. The matrix M is non-negative 
and | det(M)| = 1. One may check (cf. [6]), that in this case the corresponding Jacobian 
is given by 



J{PM n ){x) 



M K(c 1 x 1 + c 2 x 2 + ... + c a+b x a+b ) a + b ' 

where q is the sum of terms over the ith column of M and K > is a constant depending 
on the dimension a + b but independent of the matrix M. 

In order to estimate the integrals above, we have to estimate the Jacobian on the set 
0. Let x € Q. Since its coordinates are ordered in the nondecreasing order and their sum 
is equal to one, we must have x a +b > I /(a + b). Together with (6.1), we get 

1 1 



, < Xn < X a +l < . . . < X n +h < " < 1. 

2(a + b) ~ + ~ ~ a+b ~ b ~ 



This implies that on the set O we have 
1 



->a+b 



\a+b 



< J(P 



ill 



< 



(a + b) 



a+b 



K(c a + c a+1 + ... + c a+b ) a + b ' 



K(ci +c 2 + ... + c a+b ) 
We claim that 

ci + c 2 + . . . + c a -\ < (a - l)(c a + c a +i + • • . + c a+b ) 
In order to prove this assertion it is enough to show that 



(6.5) 



(6.6) 



(6.7) 



max Ci 

\<i<a+b 



max Ci 

a<i<a+b 



Recall that the matrix M is a finite product of matrices whose column vectors are 
described by (3.3). We will proceed by induction on the number of matrices involved in 
the product. If M = M V1 is just one of the elementary matrices, the claim follows directly 
from (3.3 ). Now suppose that our inequality is true for any product of elementary matrices 

>+i) 



up to the length n and let M 
over the ith column of M and c 



M n M^ 2 ---M nn M Vn+1 . Let c\ 



(») 



the sum over the ith column of M vri M 7r2 



be the sum of terms 
Once 



again from (3.3), we get that there exists a < i < a + b, such that 

> 



>+l) = » + c (n) 



o+l 



+ ... + C 



max i 

a<k<a+b 



» 



max < 

l<fc<a+fe 
(n+1) _ 



» 



» 



This easily 



Moreover, for every j ^ i there exists 1 < k < a + b such that Cj 
implies (6.7). 

Let us complete the proof of Proposition 6.1 The inequalities (6.6) and (6.7) imply 
together 

2 a+h (a + b) a+b 



1 



M 



< 



Ka a + b (c a + c a+1 + ... + c a+b ) a + b 
on 0. Thanks to this estimation we may write 

Leb(x e eS.„*^ = P n (x) £ f ) _ / f J(P M n )dx 



K(c a + c a+l + ... + c a+b ) a + b 



Leb(& { J&.. Mn 



> 



/e J(PM n )dx ~ V 2a(a + b) J Leb(Q) 



a+b 



Leb(f) 



□ 
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7 Proofs of main results 



We are now ready to prove our main results. 

Proof of Theorem \l.l\ In order to show that the set A is absorbing for the map T Q) &, it is 
enough to show that the set T is absorbing for the first return map P : — > 0. We will 
rather work in the projected space and show that the projected set T is absorbing for the 
projected first return map P. It is sufficient, since the projection map p : A a+b — > A sends 
any set of positive (a + 6)-dimensional Lebesgue measure in A a+b onto some set of positive 
(a + b — l)-dimensional Lebesgue measure in A. 

We may rewrite the decomposition (5.3) of in the projected version 



e = p-(»-i)(f)u (j j ... j e 



p>(«) 

Vl*...*V n 



ki,k2,...,k n =l gn* 1 . 



a . b 



From (6.2) we deduce that 

Leb{& \ p- n (T)) < (1 - a)Leb(Q \ p-^-^f)), 

for every n > 1. Since a G (0, 1) is a constant independent of n, this clearly implies that 
r is absorbing for P. □ 

Combining Theorem 1 1 . 1 1 and Lemma 4.5, we prove Theorem 1.2 which gives a precise 
description of the limiting behavior of the orbits. 



Proof of Theorem 1.2, Fix a > 1 and b > 2. Since the set A is absorbing for the map T a ^, 
it is enough to show that 



lim T a fe fe (x) = (0,...,0,x- 2 ,...,< 



k— >oc 



with x^j_ 2 > f° r almost every a; G A From Lemma 4.5 we know that > almost 

surely. If b = 2 we are done. 

If b > 3 we still have > for almost every x G A. Define the exceptional set 
Eb_i = {x £ A : x^ 3 +b _ 1 = 0} and let x G EV-l- It is easy to see that there exists some 
A;q > 1 such that for every k > ko we have 



(fc). 



In other terms, starting from some koth iteration, the last coordinate is the largest one 
along the whole T aj fe-orbit of x. Since x G we get 

lim x^ = lim (2tl#, • • .,*2ti),*S8 " E^ ) = (0, • • • ,0,C +6 ). 

Define E = {i£ A"^ 1 : T* 6 _ 1 (a;) ->■ (0, . . . ,0)}. Once again, from Theorem 1 1 . 1 1 and 
Lemma 4.5 we deduce that S is of zero Lebesgue measure in A a+ The analysis above 
shows that if x G C A a+b , then there exists k > 1 such that (xf o) , . . . , G E. 

We may write 

oo 

Eb-iC |J T^(Exf + ). 
fe =l 
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Since T a ^ restricted to any cylinder preserves Lebesgue measure, by Fubini Theorem we 
get that Eb_i is a set of zero measure in A a+b . Equivalently, for almost every x G A a+fe 
we have x^ +b _ 1 > 0. If b = 3 we are done. 

If b > 3 we may repeat the argument above to show that all the exceptional sets 
Eb_i = {x G A : x^ +b _ i = 0}, i = 2, . . . , b — 2, are null measure sets. □ 

It is clear from Theorem |1.2| that for every a > 1 and 6 > 2 the set D is absorbing for 
the map T a ij. Obviously, this is also the case for the projected Schweiger map S a ^ : B — y B 
defined in Section 1. 



Next we will prove Theorem 1.3 which together with Theorem 1.4 give our partial 



answer to the second Schweiger conjecture. To this end, we recall a result from [7J. 

Theorem 7.1 ([7], Corollary 8.2). The map Tn is ergodic with respect to Lebesgue mea- 
sure. 

We use it to show that Sip ■ B — > B is ergodic. 



Proof of Theorem 1.3, Recall that the transformation Si 2 is defined as a projection of the 
map Ti 2 : A 3 — > A 3 on the set B = {x G A 2 : X2 < 1}. On the projected set D = {x G 
i? : xi + X2 < 1}, which is absorbing for Sip, we have Ti^{x) = (Ti^xi, X2), £3 — x\). 
Thus i 

^i^^) = T X i(x). 

1 — %\ 

Moreover, as k — > 00, for almost every x G D we have T 1 fc 2 (x) — )• (0, 0, X3 ) where X3 > 0. 
This implies that for almost every x G D the corresponding Si^-orbit converges to the 
origin. We get that Sip is totally dissipative with respect to Lebesgue measure. 
Now consider the map (p : A 2 — > D defined by 

, xi x 2 

<P(Xl,X2) 



1 + Xl + X 2 ' 1 + Xl + X 2 

It is invertible and sends A 2 onto the set {x G D, xi + X2 < 1} of full Lebesgue measure in 
D. Moreover, c/> and c/> _1 send null measure sets to null measure sets. One may check that 



Ti a =6 1 o Si. 



o 



This together with Theorem 7.1 imply the ergodicity of Si 2 on the set D. Since D is 



almost surely absorbing, this also implies the ergodicity on the larger set B. □ 

Just as the map S12 is related to the Euclidean algorithm, for a > 2 the dynamics 
of the transformation S a p ■ D — > D is related to the Brun algorithm. Indeed, the subset 
A = {x G D : ct(x) = 1} is invariant under £^2 and its dynamics restricted to this set 
coincides with the projection of the Brun algorithm T i. Our transformation S a p may 
be thus seen as an extension of the ergodic dynamics of the Brun algorithm on A to the 
larger set {x : a(x) < 1}. It is thus natural to conjecture that its dynamics is ergodic also 
in the case a > 2. 

To show that S a b : B — > B is not ergodic for b > 3 we use the limiting behavior of 
orbits to exhibit an invariant nonconstant function. 
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Proof of Theorem 1.4 For every x £ B we have 



lim 5; 

k— >oo 



a,b\ 



i x a+b-l, 



Define / : B 



by /(*) 



It is measurable and 5 aj b-invariant. For almost every 



x E D, an argument similar to one used in the proof of Lemma 4.7 gives that 

00 _ x a+b-l ~ (Xl + ■ ■ ■ + Xg+l) 

Xa+b - 1 ~ l-( Xl + ... +Xa+1 ) ■ 

This shows that / is nonconstant. The map S a b is not ergodic for b > 3. 



□ 



8 Other subtractive algorithms 

Fix a > 2, b > 1 and 1 < % < o — 1. We may alter slightly the definition of the 
transformation T a b by subtracting the coordinate Xj instead of x a . The new transformation 
T : A a+b -)■ A a+& is given by the formula 

T(X\) X2i ■ ■ ■ 1 X a _j_f,) — 1X\X\i ■ ■ ■ 1 X ai X a -\-\ Xj, . . . , X a _|_fr Xj), 

where ir is a permutation which arranges the coordinates in ascending order. This new 
family of transformations also contains some well studied ones. For example, if i = 1 and 
b = 1 we get the Selmer algorithm. 

We may remark the following properties of the set A with respect to the transformation 
T defined above. 

Lemma 8.1. (1) For every a > 2, b > 1 and 1 < i < a — 1 we have T{A) C A. 

(2) 7/6 < o + l - i taera r -1 (A) C A. 
Proof. (1) Let x £ A and y = 7"(x). We have 

by a+b = 6max{x a ,x a+fe - Xj} > b(x a+b - Xj) > a(x) - bxt = a{y). 
This implies y G A. 
(2) Let x S c A and y = T(x). Suppose that 6 < a + 1 — i . We have 

b(x a+b - Xi) < a(x) - bxi = a(y) 

and 

a+b i—1 a a+b a+b 

bx a < V] Xj < x j + y^(. x j - x i) + x j < ^2 x j - bxi = cr{y). 

j=a+l j=l j=i j=a+l j=l 

Together this gives by a+b = &max{x a ,x Q+b - Xj} < a(y). 

□ 

In virtue of the last lemma, when looking for the globally absorbing set one has to 
assume that 

b > a + 3 — i. 

Under this assumption, we may ask if A is absorbing for the new algorithm T. The 
difficulty of this question comes from the fact that the cylinders for the map T are not 
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full as it was the case for T a ^. This is the main reason why the argument developed in 
Section 5 does not work in this case. We are no longer able to estimate the proportion of 
points absorbed by A from every cylinder. 

However, it is rather straightforward to check that the arguments presented in Section 
3, used to study limiting properties of orbits, are still valid for the new transformation T. 

Proposition 8.2. Let a > 2, b > 1 and 1 < i < a — 1. 

(1) For almost every x G A a+b the corresponding orbit converges to (xf^xEf, • • • > x ^+&) 
with xf = . . . = xf +1 = 0. 

(2) For almost every x G A we have x^ +b > 0. 

(3) For almost every x G A a+b whose orbit does not converge to the origin, there exists 
k > 1 such that T k (x) G A. 
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